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INSTITUT FRANCAIS DU PETROLE 



1 . S The applicant is hereby notified that the International Search Report has been established and is transmitted herewith. 
Filing of amendments and statement under Article 19- 

The applicant is entitled, if he so wishes, to amend the claims of the International Application (see Rule 46): 

■ ^itlrnatfonTl^yeL^ ^ '^^"^ ^^'^ transmittal of the 

international Search Report; however, for more details, see the notes on the accompanying sheet. 

Where? Directly to the International Bureau of WlPO 
34, chemin des Colombettes 
1211 Geneva 20, Switzerland 
Fascimile No.: (41-22) 740.14.35 

For more detailed Instructions, see the notes on the accompanying sheet. 
' □ »?SLrto^a^^^^^^^ ^^^^^^ ''^ --^"^^^^ '^at the declaration under 



3. □ With regard to the protest against payment of (an) additional fee(s) under Rule 40.2. the applicant is notified thaf 

□ no decision has been made yet on the protest: the applicant will be notified as soon as a decision is made. 

4. Further actlon(s): The applicant is reminded of the following: 

^ W l^^l priority date, the international application will be published bv the International Burea,, 



N ame and m ailing address of the International Searching Authority 

J European Patent Office, P.B. 5818 Patentlaan 2 
NL-2280 HV Rijswijk 
Tel. (+31-70) 340-2040, Tx. 31 651 epo nl 
Fax: (+31 -70) 340-301 6 



Authorized officer 

Jacobus Constant 
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ES TO FORM PCT/ISA/220 



Ad:"ir.".i!^^^^^^^^^^^^ '° ••^^P'ovision, of tt,e PCT, the PCT Regulations and .he PCX 



INSTRUCTIONS CONCERNING AMENDMENTS UNDER ARTICLE 19 

T>ie applicant has after having received the international search report, one opportunity to amend the claima of th« 
.ntemat.onal app(K:al.on. It should however be emphasized that, since all parts of th^ irtemiVnar^otii^^^^^^ 
tT^^"'^:^ drawings) may be amended during (he international preliminary ex^Sn p^^™ fh^^e'^'ruil^ 

^ ^?"<^'"»''«B o' 'he claims under Article 19 except where, e g, the applicant wiite tK^lTtterto ^hlth^H 

Furt*h:^oTrshou'ld L™:™h^l~Ht r "^"T '^^ ^^^^ meSr^^oTe'tnTem^ti^^L" Su^^ron"" 

funnermore, rt should be emphasized that provisional protection is available in some States only. 

What parts of the Intematlonai application may be amended? 

Under Article 1 9. only the claims may be amended. 

»h«"l^JJ^t.'"*'TD'T°'-^"'c' ''^^ <='^"«» ""^y also amended (or further amended) under Article 34 before 
%:raiSe': A^di S^"*"- '^'^ °' '"t--^*-^ ^PP--^-" -ay be amended under Article 28 



When? 



Wrthin 2 months from the dale of transmittal of the intomationaJ search r«nr.rt r.r 1 a _ *i. 

^ h;:no'^rrt"' ""d" " '^"'^ *^ ~S howr^fhaur ^ Jn^d^r^tw^^^^^^^^ 

^wS^ I °fL^ ^^7 '"".f' International Bureau after the expiration oTthe 

(^Se ^ 1) °™ '"^ completion of the technical preparations for intemat.^nlK^ion 



Where not to file the amendments? 



M^rnTo^fari^rcrrJ^*:?:]^^^^ O"*- or the 

Where a demand for international preliminary examination has been/Is filed, see below. 



^'e;;lroMl!:«d°«:ra,"^^^^^^^ °^ - — '''e text of 

^e^reTd^::* 3*.h1 s^^'SlTal^^ed''^' °' '"^^ ^'^'^^ — * °' - — *"en. or 

be renumbered consecutively (Administrative Instructions. LctioraOsTb)) renumbered, they must 

The amendments must be made In the language In which thelntematlonal application i« to be published. 

What documents must/may accompany the amendments? 
Letter (Section 205(b)): 

The amendments must be submitted with a letter. 

The letter will not be published with the intemahonal application and the amende h^i™. i> i-i . 
confused with the "Statement under Article 1 9(1)" (see'^olow?Cn^'r •SUtem%"^^ndr/^ici:^^^^^^^^^^ 
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^^NOTES TO FORM PCT/iSA/220 (coni^|ij) 



Th© letter must indicate the differences between the claims as filed and the claims as amended. It must, in 
particular, indicate, in connection with each claim appearing in the international application (it being underatood 
that identical indications concerning several claims may be grouped), whether 

(i) the claim is unchanged; 

(li) the claim is cancelled; 

(iii) the claim is new; 

(iv) th© claim replaces one or more daims as filed; 

(v) th© claim is the result of the division of a claim a3 filed. 



The followtng examples Illustrate the manner in which amendments must be explained In the 
accompanying letter: 



1 [Where onginally there were 48 claims and after amendment of some claims there are 51 ]• 

■aajms 1 to 29, 31 , 32. 34, 35, 37 to 48 replaced by amended claims beanng the same numbers- 
daims 30, 33 and 36 unchanged; new claims 49 to 51 added." 

2. [Where originally there were 1 5 claims and after amendment of all daims there are 11 ]: 
"Claims 1 to 15 replaced by amended claims 1 to 1 1 .* 

3. [Where originally there were 14 claims and the amendments consist in cancelling some claima and in adding 
new claims]; ^ 
•Claims 1 to 6 and 1 4 unchanged; claims 7 to 13 cancelled; new claims 15, 16 and 17 added." or 

"Claims 7 to 13 cancelled; new claims 15, 16 and 17 added; all other claims unchanged." 

4. [Where various kinds of amendments are made): 

■Claims 1-10 unchanged; claims 1 1 to 1 3, 1 8 and 1 9 cancelled; claims 1 4, 1 5 and 1 6 replaced by amended 
daim 14; claim 1 7 sutsdivided into amended claims 15, 16 and 17; new daims 20 and 21 added." 



"Statement under article 19(1)" (Rule 46.4) 

The amendments may be accompanied by a statement explaining th© amendments and indicating any impact 
that such amendments might have on the description and the drawings (which cannot be amended under 
Artid© 19(1)). 

The statement will be published with the international application and the amended claims. 
It must be In the language In which the International appplicatton is to be published. 

It must be brief, not exceeding 500 words if in English or if translated into English. 

It should not be confused with and does not replace the letter indicating the differences between the claims 
as filed and as amended, tt must be filed on a separate sheet and must be identified as such by a heading 
preferably by using the words "Statement under Article 1 9(1)." 

tt may not contain any disparaging comments on th© international search report or the relevance of citations 
contained in that report. Reference to citations, relevant to a given daim. contained in the intemationaJ search 
report may be made only in connection with an amendment of that claim. 



Consequence if a demand for international preliminary examination has already been filed 

If, at the time of filing any amendments under Article 1 9, a demand for international preliminary examination 
has already been submitted, the applicant must preferably, at the same time of filing the amendments with the 
International Bureau, also file a oopy of such amendments with the International Preliminary Examining 
Authority (see Rule 62.2(a). first sentence). 



Consequence with regard to translation of the tntematlonal application for entry Into the national phase 

TTie applicant's attention is drawn to the fact that, where upon entry into the national phase, a translation of the 
daims as amended under Article 1 9 may have to be furnished to the designated/elected Offices, instead of or 
in addition to, the translation of the daims as filed. 

^further details on the requirements of each designated/elected Office, see Volume H of the PCT Applicant's 
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ATENT COOPERATION TRE 



PCT 



INTERNATIONAL SEARCH REPORT 

(PCT Article 1 8 and Rules 43 and 44) 



Applicant's or agent's file reference 

4530/JC 


pQp FURTHER ®^ Notification of Transmittal of International Search Report 

(Form PCT/ISA/220) as well as, where applicable, item 5 below. 

ACTION 


International application No. 

PCT/FR 00/01853 


International filing date (day/month/year) 

30/06/2000 


(Earliest) Priority Date (day/month/year) 

02/07/1999 


Applicant 

INSTITUT FRANCAIS DU PETROLE 



This International Search Report has been prepared by this International Searching Authority and is transmitted to the applicant 
according to Article 18. A copy is being transmitted to the International Bureau. 

This International Search Report consists of a total of 3 sheets. 

It is also accompanied by a copy of each prior art document cited in this report. 



1 . Basis of the report 

a. With regard to the language, the international search was carried out on the basis of the international application in the 
language in which it was filed, unless otherwise indicated under this item. 

I I the international search was carried out on the basis of a translation of the international application furnished to this 
Authority (Rule 23.1(b)). 

b. With regard to any nucleotide and/or amino acid sequence disclosed in the international application, the international search 
was carried out on the basis of the sequence listing ; 
I I contained in the international application in written form. 

filed together with the international application in computer readable form, 
furnished subsequently to this Authority in written form, 
furnished subsequently to this Authority in computer readble form. 



2. 
3. 



□ 
□ 
□ 
□ 

□ 



□ 
□ 



the statement that the subsequently furnished written sequence listing does not go beyond the disclosure in the 
international application as filed has been furnished. 

the statement that the information recorded in computer readable form is identical to the written sequence listing has been 
furnished 

Certain claims were found unsearchable (See Box 1). 
Unity of Invention Is lacking (see Box II). 



4. With regard to the title, 

pT] the text is approved as submitted by the applicant. 

I I the text has been established by this Authority to read as follows: 



With regard to the abstract, 

[X] the text is approved as submitted by the applicant, 
the text has been established, according to Rule 3J 

within one month from the date of mailing of this international search report, submit comments to this Authority. 
The figure of the drawings to be published with the abstract is Figure No. 1 



j — I the text has been established, according to Rule 38.2(b), by this Authority as it appears in Box III. The applicant may, 



PC] as suggested by the applicant. Q None of the figures. 

I I because the applicant failed to suggest a figure. 

I I because this figure better characterizes the invention. 
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INTERNATIONAL SEARCH REPORT 



International Application No 

T/FR 00/01853 



A. CLASSIFICATION OF SUBJECT MATTER 
IPC 7 GOlVll/OO 



According to Internationa! Patent Classification (IPC) or to both national classification and IPC 



B. FIELDS SEARCHED 



Minimum documentation searched (ciasstfication system followed by classification symbols) 

IPC 7 GOIV 



Documentation searched other than minimum documentation to the extent that such documents are included in the fields searched 



Electronic data base consulted during the international search (name of data base and, where practical, search terms used) 

WPI Data, INSPEC, COMPENDEX, EPO-Internal 



0. DOCUMENTS CONSIDERED TO BE RELEVANT 



Category ° Citation of document, with indication, where appropriate, of the relevant passages 



Relevant to claim No. 



EP 0 911 652 A (WESTERN ATLAS INT INC) 
28 April 1999 (1999-04-28) 

column 2, line 57 -column 3, line 13 

column 4, line 34 -column 7, line 12 

US 5 838 634 A (HELWICK JR STERLING J 
AL) 17 November 1998 (1998-11-17) 
column 19, line 30 -column 21, line 30 



ET 



DEUTSCH C V ET AL: "GEOSTATISTICAL 
TECHNIQUES IMPROVE RESERVOIR MANAGEMENT" 
PETROLEUM ENGINEER INTERNATIONAL, US, HART 
PUBLICATIONS, 

vol. 69, no. 3, 1 March 1996 (1996-03-01), 
pages 21-22,24-27, XP000596614 
ISSN: 0164-8322 

page 25, left-hand column -page 26, middle 
column 

-/-- 



1-4 



1-4 



1-4 



m 



Further documents are listed in the continuation of box C. 



|)( j Patent family members are listed in annex. 



** Special categories of cited documents : 

"A" document defining the general state of the art which is not 
considered to be of particular relevance 

"E" earlier document but published on or after the international 
filing date 

"L" document which may throw doubts on priority ctaim(s) or 
which is cited to establish the publication date of another 
citation or other special reason (as specified) 

"O" document referring to an oral disclosure, use, exhibition or 
other means 

"P" document published prior to the international filing date but 
later than the priority date claimed 



"T" later document published after the international filing date 
or priority date and not in conflict with the application but 
cited to understand the principle or theory underlying the 
invention 

"X' document of particular relevance; the claimed invention 
cannot be considered novel or cannot be considered to 
involve an inventive step when the document is taken alone 

'Y" document of particular relevance; the claimed invention 

cannot be considered to Involve an inventive step when the 
document is combined with one or more other such docu- 
ments, such combination being obvious to a person skilled 
in the art. 

■&■ document member of the same patent family 



Date of the actual completion of the international search 



6 September 2000 



Date of mailing of the intemational search report 



13/09/2000 



Name and mailing address of the ISA 

European Patent Office. P.B. 5818 Patentiaan 2 
NL - 2280 HV Riiswiik 
Tel. (+01-70) 340-2040. Tx. 31 651 epo nl. 
Fax: (+01-70)340-^16 



Authorized officer 



Swart jes, H 
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INTERNATIONAL SEARCH REPORT 



^ynter 



iternatlonal Application No 

T/FR 00/01853 



C.(Continuatlon) DOCUMEI^S CONSIDERED TO BE RELEVANT 



Category * 



Citation of document, with Indication, where appropriate, of the relevant passages 



Relevant to claim No. 



CARR JAMES R ET AL: "Use of geostati sti cs 
for accurate mapping of earthquake ground 
motion " 

GEOPHYSICAL JOURNAL, US , GORDON AND BREACH, 
NEW YORK, NY, 

vol. 97, 1 January 1989 (1989-01-01), 
pages 31-40, XP002097135 
ISSN: 0275-9128 



FofTTi PCT/ISA/210 {continuation of second sheet) (July 1992) 
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INTERNATIONAL SEARCH REPORT 

nation on patent family members 



Patent document 
cited in search report 



Publication 
date 



Uiternational Application No 

CT/FR 00/01853 



Patent family 
member(s) 



Publication 
date 



EP 0911652 



28-04-1999 



US 
NO 



5995906 A 
984499 A 



30-11-1999 
06-04-1999 



US 5838634 ' 


A 


17-11-1998 


CA 


2251365 


A 


16-10-1997 








EP 


0891562 


A 


20-01-1999 








NO 


984621 


A 


04-12-1998 








WO 


9738330 


A 


16-10-1997 
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PC 



REQUETE 



Le soussigne requiert que la presente demande 
intemationale soit traitee conformement au Traite de 
cooperation en matiere de brevets. 



ve a r office recepteur 



Demande Internationale n** 



Date du depot international 



Nom de Toffice recepteur et "Demande internationale PCT" 



Reference du dossier du depos^int ou du mandataire (facultatij) 
(12 caracteres au maximum) 4530 / JC 



Cadre n'' I TITRE DE L'INVENTION METHODE POUR DEFORMER GRADUELLEMENT DES SIMULATIONS 
SEQUENTELLES D'UN MILIEU HETEROGENE TEL QU'UNE ZONE SOUTERRAINE 



Cadre II DEPOSANT 



Nom et adresse : (Nom de famille suivi du prenom; pour une personne morale, designation 
officielle complete. L 'adresse doit comprendre le code postal et le nom du pays. Le pays de 
I adresse indiquee dans ce cadre estl'Etat oil le deposant a son domicile si ducun domicile 
n 'est indique ci-dessous.) 

INSTITUT FRANgAIS DU PETROLE 
1 & 4 avenue de Bois Preau 
92852 Rueil Malmaison cedex 
FRANCE 



□ 



Cette personne est aussi 
inventeur. 



n" de telephone 

01 47 52 60 00 



n" de telecopieur 

01 47 52 70 03 



n" de teleimprimeur 



Nationalite (nom de TEtat) : FRANCE 


Domicile (nom de I'Etat) : FRANCE 


Cette personne est i 1 tousles Etats 

deposant pour : | | designes 


"vTI tous.Ies Etats designes sauf 1 1 les Etats-Unis d'Amerique | 1 les Etats indiques dans 

•^1 les Etats-Unis d'Amerique | | seulement | | le cadre supplementaire 


Cadre n° III AUTRE(S) DEPOSANT(S) OU (AUTRE(S)) INVENTEUR(S) 



Nom et adresse : (Nom de famille suivi du prenom; pour une personne morale, designation 
officielle complete. L 'adresse doit comprendre le code postal et le nom du pays. Le pays de 
I adresse indiquee dans ce cadre est I'Etat oil le deposant a son domicile si ducun domicile 
n 'est indique ci-dessous.) 

ELF EP 
Tour Elf 

2, place de la Coupole, La Defense 6 

92400 Courbevoie 

FRANCE 


Cette personne est ; 

deposant seulement 

[ 1 deposant et inventeur 

1 1 inventeur seulement 
(Si cette case est cochee, 
ne pas remplir la suite.) 


Nationalite (nom de I'Etat) : 

FRANCE 


Domicile (nom de TEtat) : 

FRANCE 



Cette personne est 
deposant pour : 



□ tous les Etats 
designes 



tou^ les Etats designes sauf 
les Etats-Unis d'Amerique 



nles Etats-Unis d'Amerique I 1 les Etats indigues dans 
seulement | | le cadre supplementaire 



I I D'autres deposants ou inventeurs sont indiques sur une feuille annexe. 



Cadre n"* IV MANDATAIRE OU REPRESENT ANT COMMUN; OU ADRESSE POUR LA CORRESPONDANCE 



La personne dent I'identite est donnee ci-dessous est/a ete designee pour agir au nom du ou | i mandataire fvl representant commun 
des deposants aupres des autorites mternationales competentes, comme: 1 | U^l ^ 


Nom et adresse : (Nom de famille suivi du prenom; pour une personne morale, designation officielle 
complete. L 'adresse doit comprendre le code postal et le nom du pays.) 

ELMALEH Alfred 

INSTITUT FRANCAIS DU PETROLE 

1 & 4, Avenue de Bois Preau 

92852 RUEIL MALMAISON CEDEX FRANCE 


n* de telephone 

01 47 =^9 ft? 79 


n* de telecopieur 

01 47 52 70 03 


n" de teleimprimeur 



□ 



Adresse pour la correspondance : cocher cette case lorsque aucun mandataire ni representant commun n'est/n*a ete designe 
et que Tespace ci-dessus est utilise pour indiquer une adresse speciale a laquelle la correspondance doit etre envoyee. 



Formulaire PCT/RO/101 (premiere feuille) (juiHet 1998; reimpression Janvier 2000) Voir les notes relatives au formulaire de requete 
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Feuilie n'* 



Suite du cadre n" III AUTRE(S) DEPOSANT(S) OU (AUTRE(S)) INVENTEUR(S) 



Si aucun des sous-cadres suivants n 'est utilise, cette feuilie ne doit pas etre incluse dans la requete. 


Nom et adresse : (Nam de famille siiivi du prenom; pour une personne morale, designation 
offtcielle compleie. L 'adresse doU comprendre le code postal ei le nom du pays. Le pays de 
I adresse indiquee dans ce cadre est I 'Etat oii le deposant a son domicile si aucun domicile 
n 'est indique ci-dessous.) 

HU Lin-Ying 

4 rue Louis Borthou 

64140 Billere 

FRANCE 


Cette personne est : 

I [ deposant seulement 

deposant et inventeur 

I 1 inventeur seulement 
(Si cette case est cochee, 
ne pas remplir la suite.) 


Nationalite (nom de TEtat) : FRANCE 


Domicile (nom de I'Etat) : FRANCE 


Cette personne est 1 1 tous lesEtats 1 | tous les Etats designes sauf [wf les Etats-Unis d'Amerique 1 1 les Etats indiques dans 

deposant pour : I 1 designes 1 1 les Etats-Unis d'Amerique seulement 1 1 le cadre supplementaire 


Nom et adresse : (Nom de famille suivi du prenom; pour une personne morale, designation 
offtcielle complete. L 'adresse doit comprendre le code postal et le nom du pays. Le pays de 
I adresse indiquee dans ce cadre est I 'Etat oil le deposant a son domicile si aucun domicile 
n 'est indique ci-dessous.) 

NOETINGER Benoit 
57 rue Emile Guichenne 
64000 Pau 
FRANCE 


\ 

Cette personne est : 

1 I deposant seulement 

[^>^ deposant et inventeur 

1 1 inventeur seulement 
(Si cette case est cochee, 
ne pas remplir la suite.) 


Nationalite (nom de TEtat) : FRANCE 


Domicile (nom de TEtat) : FRANCE 


Cette nersonne est 1 1 tous les Etats 1 1 tous, les Etats designes sauf les Etats-Ums d'Amerique 1 1 les Etats indiques dans 

de^posant p^^^^^ 1 1 designes \ | les Etats-Unis d'/Siierique X seulement | 1 le cadre suppfementaire 


Nom et adresse : (Nom de famille suivi du prenom: pour une personne morale, designation 
officielle complete. L 'adresse doit comprendre le code postal et le nom du pays. Le pays de 
I adresse indiquee dans ce cadre est I 'Etat oil le deposant a son domicile si aucun domicile 
n 'est indique ci-dessous.) 


Cette personne est : 

I 1 deposant seulement 

I 1 deposant et inventeur 

1 1 inventeur seulement 
(Si cette case est cochee. 
ne pas remplir la suite.) 


Nationalite (nom de TEtat) ; 


Domicile (nom de l*Etat) : 


Cette personne est 1 1 tous les Etats 1 1 tous les Etats designes sauf 1 1 les Etats-Unis d'Amerique 1 | les Etats indiques dans 

deposant pour 1 1 designes | | les Etats-Unis d'/5nerique | | seulement | 1 le cadre suppfementaire 


Nom et adresse : (Nom de famille suivi du prenom; pour une personne morale, designation 
officielle complete. L 'adresse doit comprendre le code postal et le nom du pays. Lepays de 
I adresse indiquee dans ce cadre est I'Etat oil le deposant a son domicile si aucun domicile 
n 'est indique ci-dessous.) 


Cette personne est : 

I 1 deposant seulement 

1 1 deposant et inventeur 

1 1 inventeur seulement 
(Si cette case est cochee, 
ne pas remplir la suite.) 


Nationalite (nom de I'Etat) : 


Domicile (nom de I'Etat) : 


Cette personne est 1 1 tous les Etats | 1 tous. les Etats designes sauf 1 1 les Etats-Unis d'Amerique 1 1 les Etats indi<jues dans 

de>)sant pou^ : 1 1 designes | | les Etats-Unis d'.5nerique | 1 seulement | 1 ie cadre supplementaire 


[ 1 D'autres deposants ou inventeurs sont indiques sur une autre feuilie annexe. 
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Cadre n° V DESIGNATION 



Feuille n** 



'ATS 



L^s desianations suivantes sont faites conformement a la regie 4. 9. a) (cocher les cases appropriees; une au mains doitl'etre) : 
Brevet regional 

□ AP Brevet ARIPO : GH Ghana. GM Gamble, KE Kenya, LS Lesotho, MW Malawi, SD Soudan, SL Sierra Leone, 

SZ Swaziland, TZ Republique-LFnie de Tanzanie, UG Oug'anda, ZW Zimbabwe et tout autre Etat qui est un Etat contractant du 
Protocole de Harare et du PCX 

□ EA Brevet eurasien : AM Armenie,AZ Azerbaidjan, BY Bdarus, KG Kirghi2istan,KZ Kazakhstan, MDRepubliquede Moldova, 

RU Federation de Russie,TJ Tadjikistan, TM Turkmenistan et tout autre Etat qui est un Etat contractant de la Convention sur 
le brevet eurasien et du PCT 

S EP Brevet europeen : AT Autriche, BE Belgique, CH et LI Suisse et Liechtenstein, CY Chypre, DE AUemagne, 
DK Danemark, ES Espasne, FI Finlande, FR France, GB Royaume-Uni, GR Grece, IE Irlande, IT Italie, 
LU Luxembourg, MC Monaco, NL Pays-Bas, PT Portugal, SE Suede et tout autre Etat qui est un Etat contractant de la 
Convention sur le brevet europeen et du PCT 

□ OA Brevet OAPI : BF Burkina Faso, BJ Benin, CF Republique centrafricaine, CG Congo, CI Cote d'lvoire, 

CMCameroun, GA Gabon, GN Guinee, GW Guinee-Bissau, ML Mali, MR Mauritanie, NE Niger, SN Senegal, 
TD Tchad, TG Togo et tout autre Etat qui est un Etat membre de TOAPI et un Etat contractant du PCT (si une autre forme 

de protection ou de traitement est souhaiiee, le preciser sur la ligne pointille^ 

Brevet national (si une autre forme de protection ou de traitement est souhaitee, le preciser sur la ligne pointillee) : 

□ AE Emirats arabes unis □ LR Liberia 

□ AL Albanie □ 

□ AM Armenie C] 

□ AT Autriche □ 

Austral ie D 



□ AU 

□ AZ 

□ BA 

□ BB 

□ BG 

□ BR 

□ BY 

□ CA 



Azerbaidjan 



LS 
LT 
LU 
LV 
MA 



Lesotho . . 
Lituanie 
Luxembourg 
Lettonie 
Maroc . . . . 



Bosnie-Herzegovine □ MD Republique de Moldova 



Barbade 

Bulgarie 

Bresil 

Belarus □ 

Canada □ 

□ CH et LI Suisse et Liechtenstein D 

□ CN Chine IS 

Costa Rica CD 

Cuba □ 



[U MG Madagascar 

□ MK Ex-Republique yougoslave de Macedoine 



□ CR 

□ CU 

□ CZ 

□ DE 

□ DK 
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TZ 
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Ouzbekistan . . . 
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Afrique du Sud 



Kenya □ 

Kirghizistan □ 

Republique populaire democratique de Coree . CD 

□ ZW Zimbabwe 

□ KR Republique de Coree Cases reservees pour la designation d'Etats qui sont devenus parties 

□ KZ Kazakhstan apres la publicatiorTde la presente feuille : 

□ LC Sainte-Lucie □ 

□ LK Sri Lanka □ 

Declaration concernant les designations de precaution : outre les designations faites ci-dessus, le deposant fait aussi conformement 
a la regie 4.9. b) toutes les designations qui seraient autorisees en vertu du PCT, a Texception de toute designation indiqueedans le cadre 
suppllmentaire comme etant exclue de la portee de cette declaration. Le deposant declare que ces designations additionnelles sont 
faites sous reserve de confirmation et que toute designation qui n'est pas confirmee avant Texpiration d'un delai de 15 mois acompter 
de la date de priorite doitetre consideree comme retiree par le deposant a Texpirationde ce delai. (La confirmation (y compris les taxes) 
doit parvenir a I 'office recepteur dans le delai de 15 mois.) 
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Cadre n° VI 



REVENDICA 



Feuille n** 



DE PRIORITE 



D'autres revendications de priori te sont 
indiquees dans le cadre supplementaire. 



Date de depot 
de la demande anterieure 
(jour/mois/annee) 


IN Ulllwl U 

de la demande anterieure 


Lorsque la demande anterieure est une : 


demande nationale : 
pays 


demande regionale :* 
office regional 


demande internationale : 
office recepteur 


(1) 

OZjuillet 1999 
m9/n7/QQ^ 


99/08.605 


FRANCE 
















(3) 












I presente 



* Si la demande anterieure est une demande ARIPO, il est obligatoire d'indiquer dans le cadre supplementaire au mains un pays pa-tie a la Convention 
de Paris pour la protection de la propriete industnelle pour lequel cette demande anterieure a ete deposee (regie 4.10.b)ii}). Voir le cadre supplementaire. 

Cadre n° VII ADMINISTRATION CHARGEE DE LA RECHERCHE INTERNATIONALE 



Choix de T administration chargee de la recherche 
internationale (ISA) (si plusieurs administrations 
chargees de la recherche internationale sont competentes 
pour proceder a la recherche internationale, indiquer 
I 'administration choisie; le code a deux lettres peut etre 
utilise) : 

ISA/ 



Demande d'utilisation des resultats d'une recherche anterieure; mention de 
cette recherche (si une recherche anterieure a ete effectuee par i administration 
chargee de la recherche internationale ou demartdee a cette demiere) 

Date (jour/mois/annee) Numero Pays (ou office regional) 

12 Mai 2000 FA 580327 FRANCE 



Cadre n° VIII BORDEREAU; LANGUE DE DEPOT 



La presente demande internationale contient 
le nombre de feuilles suivant : 



requete 

description (sauf partie reservee 
au listage des sequences) 

revendications 

abrege 

dessins 

partie de la description reservee 
au listage des sequences : 

Nombre total de feuilles 



13 
2 
1 
3 
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Le ou les elements coches ci-apres sont joints a la presente demande internationale : 

1 . □ feuille de calcul des taxes 

2. □ pouvoir distinct signe 

3 . □ copie du pouvoir general; numero de reference, le cas echeant : 

4. □ explication de 1' absence d'une signature 

5. □ document(s) de priorite indique(s) dans le cadre n** VI au(x) point(s) : 

6. □ traduction de la demande internationale en (langue) : 

7. □ indications separees concernant des micro-organismes ou autre materiel 

biologique deposes 

8. □ listage des sequences de nucleotides ou d'acides amines sous forme 

dechiffrable par ordinateur 

9. B autrcs elements (preciser) : Cop^ ^Cyi^^cV (j^. reiinQTcVlg 



Figure des dessins qui 

doit accompagner 1' abrege : <X 


Langue de depot de la 
demande internationale : 


Cadre n° IX SIGNATURE DU DEPOSANT OU DU MANDATAIRE 



INSTTTUT FRANQAIS DU PETROLE 



1 et 4 avenue de Bois-Preau 
92852 RUEIL-MALMAISON CEDEX 
(France) 



INSTITUT FRAMCAIS DU PETROLE 




Alfred ELMALEH 
Chef du Departement Brevets 



1 . Date effective de reception des pieces supposees 
constituer la demande internationale : 


2. Dessins : 
1 1 re^us : 

1 1 non re9us : 


3. Date effective de reception, rectifiee en raison de la reception ulte- 
rieure, mais dans les delais, de documents ou de dessins completant ce 
qui est suppose constituer la demande internationale : 


4. Date de reception, dans les delais, des corrections 
demandees selon Particle 1 1.2) du PCT : 


5. Administration chargee de la recherche a / 
internationale (si plusieurs sont competentes) : 1 


6. 1 1 Transmission de la copie de recherche differee 
1 1 jusqu'au paiement de la taxe de recherche. 



Date de reception de Texemplaire 
original par le Bureau international : 



Reserve au Bureau international 
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FEUILLE DE CALCUL DES TAXES 
Annexe tie la requete 



Reference du dossier du 

deposant ou du mandataire AFF. 4530/JC 



Reserve a I'office recepteur 



Demande internationaien^ 



Timbre a date de Toffice recepteur 



Deposant 



INSTITUT FRANQAIS DU PETROLE et ELF EP 



CALCUL DES TAXES PRESCIUTES 

1. TAXE DE TRANSMISSION . . . 



400,0C T 



TAXE DE RECHERCHE 

Recherche internationale a effectuer par 

(Si plusieurs administrations chargees de la recherche internationale sont 
competentes en ce qui concerne la demande internationale, inscrire l& nom de celle 
qui est choisie pour la recherche internationale.) 



619879|T 



3. TAXE INTERNATIONALE 
Taxe de base 

La demande internationale contient 
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feuilles. 



30 premieres feuilles |_ 



F. 2682,86 



bl 



feuilles suivantes 



montant additionnel 



b2 



Additionner les montants portes dans les cadres 
bl et b2 et inscrire le total dans le cadre B . . . 



Taxes de designation 

La demande Internationale contient 



designations. 



F. 2682,86 



B 



577.24 



montant de la taxe de designation 



F. 1731,72 



nombre de taxes de 
designation dues (maximum 8) 

Additionner les montants portes dans les cadres B et D, et 

inscrire le total dans le cadre I 

(Les deposants de certains E tats ont droit a une reduction de 75 % sur la taxe 
internationale. Lorsque le deposant a (ou tous les deposants ont) droit a cette 
reduction, la somme devant jigurer sous I est egale d 25 % de la somme des 
montants figurant sous B et DJ 

4. TAXE AFFERENTE AU DOCUMENT DE PRIORITE (le cas echeant) . 

5. TOTAL DES TAXES DUES 

Additionner les montants portes dans les cadres 

T, S, I et P, et inscrire le resultat dans le cadre TOTAL 



4414,56 I 



100, OQ 



F. 11113,37 



TOTAL 



[ I Les taxes de designation seront payees ulterieurement. 



MODE DE PAIEMENT 

I \/ autorisation de debiter un compte I I traite bancaire 

' — ^ de depot (voir ci-dessous) ' 

I I cheque Q especes 

[ I mandat postal | | timbres fiscaux 



I I coupons 

I I autres (preciser): 



AUTORISATION CONCERNANT UN COIVIPTE DE DEPOT ^fes offices recepteurs ne permettent pas tous I'utilisation de ce mode de paiement) 
L* office recepteur/ 



[ )^ est autorise a debiter mon compte de depot du total des taxes indique ci-dessus. 

□ 
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□ 



(cette case ne peut etre cochee que si les conditions relatives aux comptes de depot etablies par I office 
recepteur le permettent) est autorise a debiter mon compte de depot de tout montant manquant - ou a le 
crediter de tout excedent - dans le paiement du total des taxes indique ci-dessus. 

est autorise a debiter mon compte de depot du montant de la taxe afferente a Tetablisg 
de priorite et a sa transmission au Bureau international de TOMPL 

mercredi 28 juin 2000 Alfred ELMALEH, 



ent du document 



Numero du compte de depot 



Date (jour/mois/annee) 



Signature 
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INSTITUT FRANgAIS DU PETROLE 
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R6f6rence du dossier du deposant ou du mandataire 
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AVIS IMPORTANT 


Demande intemationale no 
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Date du dSpdt international Oour/mois/annee) 
30 juin 2000 (30.06.00) 


Date de priority Uour/mois/annee) 
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INSTITUT FRANQAIS DU PETROLE etc 
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2. Les offices deslgnes suivants ont renonc6 d I'exigence selon laquelle cette communication doit fetre effectu6e a cette date: 

EP,NO 

La communication sera effectu6e seulement sur demande de ces °«ices De f^us. le d6posant n'est pas tenu de remettre 
de copie de la demande intemationale aux offices en question (regie 49.1)a-bis)). 
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compter de la date de pnorit6. 

II appartient exclusivement au deposant de veiller au respect du delai de 19 mois. 

II est a noter que seul un deposant qui est ressortissant d'un Etat contractant du PCT Ii6 par le chap.tre II ou qu, y a son 
domicile peut presenter une demande d'examen pr6liminaire international. 

RAPPEL CONCERNANT UOUVERTURE DE LA PHASE NATIONALE (article 22 ou 39.1)) 

de chaque office designe ou 6lu. 
du Guide du deposant du PCT. 
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121 1 Geneve 20, Suisse 
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(54) Title: METHOD FOR GRADUALLY DEFORMING SEQUENTIAL SIMULATIONS OF A HETEROGENEOUS ENVI- 
RONMENT SUCH AS AN UNDERGROUND ZONE 

(54) Titre: METHODE POUR DEFORMER GRADUELLEMENT DES SIMULATIONS SEQUENTIELLES D'UN MILIEU HE- 
TEROGENE TEL QU'UNE ZONE SOUTERRAINE 



QO 



O 



(57) Abstract: The invention concerns a method for gradually 
defonning representations or productions, generated by sequen- 
tial simulation, of a stochastic model not necessarily Gaussian of a 
physical quantity 2 in a meshed heterogeneous environment, so as 
to adjust them to a set of data concerning the structure or the state 
of the environment collected by prior measurements and observa- 
tions. It consists essentially in applying an algorithm gradually 
deforming a stochastic model to a Gaussian vector with N mu- 
tually independent variables, which is linked to a uniform vector 
with N mutually independent variables by the function of Gauss- 
ian distribution so as to define productions of the uniform vec- 
tor, and in using said productions to generate representations of 
said physical quantity z, which is set relative to the data. The in- 
vention is useful for example to display the statistic configuration 
of a quantity: the permeability of an undergroimd deposit, atmo- 
spheric pollution and others. 

(57) Abrege: Meihode pour deformer graduellement les 
repr^entations ou realisations, g^nerees par simulation 
sequentielle, d*un modele stochastique non necessairement 
gaussien d'une grandeur physique z dans un milieu heterogene 
mailld, afin de les ajuster a un ensemble de donnees relatives a la structure ou Tetat du milieu qui sont collectees par des mesures 
et observations prealables. Elle comporte essentiellement TappUcation d'un algoritiune de deformation graduelle d'un modele 
stochastique a un vecteur gaussien a N variables mumellement independantes, qui est leUe a un vecteur uniforme a N variables 
uniformes mumellement independantes par la fonction de repartition gaussieime de fa?on a d^finir des realisations du vecteur 
uniforme, et r utilisation de ces realisations pour generer des representations de cetie grandeur physique z, que Ton cale par rappon 
aux doimees. Applications par exemple pour visualiser la configuration statistique d'une grandeun la permeabilite d'un gisement 
souterrain, la pollution atmosph6rique, etc. 
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En ce qui concerne les codes a deux letires et autres abrevia- 
tions. se referer aux "Notes explicatives relatives aux codes et 
abreviations" figurant au debut de chaque numero ordinaire de 
la Gazette du PCX 
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METHODE POUR DEFORMER GRADUELLEMENT DES 
SIMULATIONS SEQUENTIELLES D'UN MILIEU HETEROGENE 
TEL QU'UNE ZONE SOUTERRAINE 

5 La presente invention a pour objet une methode pour deformer 

graduellement des representations ou realisations, generes par simulation 
sequentielle, d'un modele stochastique non necessairement gaussien d'un 
milieu heterogene, en se basant sur un algorithme de deformation graduelle de 
modeles stochastiques gaussiens. 

10 La methode selon I'invention trouve des applications dans la 

modelisation de zones souterraines ou il s'agit de generer des representations 
montrant comment est distribute une certaine grandeur physique dans une 
zone du sous-sol (la permeabilite z par exemple), qui soient compatibles au 
mieux avec des donnees observees ou mesurees : donnees geologiques, 
15 enregistrements sismiques, mesures obtenues dans des puits notamment 
mesures des variations au cours du temps de la pression et du debit de fluides 
issus d'un gisement etc. 

Etat de la technique 

Par la demande de brevet FR 98/09.018, on connait une methode pour 
deformer graduellement un modele stochastique (de type gaussien ou 
apparente) d'un milieu heterogene tel qu'une zone souterraine, contraint par 
une ensemble de parametres relatifs a la structure du milieu. Cette methode 
comporte le tirage d'un nombre p (p=2 par exemple) de realisations (ou 
representations) independantes du modele ou d'une partie au moins du modele 
choisi du milieu, parmi T ensemble de toutes les reahsations possibles et une 



20 



25 
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ou plusieurs etapes iteratives de deformation graduelle du modele en 
effectuint une ou plusieurs combinaisons lineaires successives de p 
realisations initiales independantes entre elles puis des realisations composites 
successivement obtenues eventuellement avec de nouveaux tirages, etc., les 
5 coefficients de cette combinaison etemt tels que la somme de leurs carres est 
egale a 1 . 

Les modeles gaussiens ou apparentes conviennent pour modeliser des 
champs de grandeurs continues et, de ce fait, sont mal adaptes a la 
modelisation de zones traversees par des reseaux de fractures ou des systemes 
10 de chenaux, par exemple. 

On connait par ailleurs les algorithmes dits de simulation sequentielle, 
les plus largement utilises parmi les algorithmes de simulation geostatistique. 
Bien qu'ils conviennent particulierement pour la simulation des modeles 
gaussiens, il n'impliquent pas en principe une limitation a ce type de modele. 

15 On forme une representation geostatistique d'une zone souterraine par 

exemple en la subdivisant par un reseau a N mailles et en determinant un 
vecteur aleatoire a N dimensions Z = (Zi, Z2,.... Zn) correspondant au mieux a 
des mesures ou observations obtenues sur la zone. Comme le montre par 
exemple Johnson, M.E. ; in « Multivariate Statistical Simulation » ; Wiley & 

20 Sons, New York, 1987, cette approche reduit le probleme de la creation d'un 
vecteur a N dimensions a une serie de N problemes mono-dimensionnels. Un 
tel vecteur aleatoire n'est ni necessairement multi-gaussien ni stationnaire. La 
simulation sequentielle de Z implique d'abord la definition d'un ordre selon 
lequel les N elements (Zi, Z2,.... Zn) du vecteur Z sont generes Tun apres 

25 r autre. Sans rien perdre de generality, on suppose que les N elements de Z 
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sont generes en sequence depuis Zj jusqu'aZN. Pour tirer une valeur de chaque 
element Zi, (i = 1, N), il faut realiser les operations suivantes : 

a) construire la distribution de Zj conditionnee par (Zj, Z2.. Zj.i) 
F,U,) = P(Z, <z, /Z,,Z3,....Z,_J ;et 

b) tirer une valeur de Zj de la distribution Fc(Zi). 

Dans la pratique de la geostatistique, la simulation sequentielle est 
utilisee frequemment pour engendrer des vecteurs multi-gaussiens et des 
vecteurs d'indicateur non-gaussiens. La tache principale de la simulation 
sequentielle est la determination des distributions conditionnelles Fc(Zi) (i = 1, 
N). Des algorithmes et logiciels pour esdmer ces distributions sont decrits 
par exemple dans : 

- Deutsch, C.V. et al, « GSLIB (Geostatistical Software Library) and User's 
Guide » ; Oxford University Press, New York, Oxford 1992. 

En ce qui conceme le tirage des valeurs a partir de la distribution Fc(Zi), 
il existe aussi un vaste ensemble d' algorithmes connus. 

On considere la methode de distribution inverse par laquelle on obtient 
une realisation de Zi : z- = F~^(u^), oxx est tire d'une distribution uniforme 
entre 0 et L II s'ensuit qu'une realisation du vecteur Z correspond a une 
realisation du vecteur U dont les elements Uj, U2, ....,Un, sont mutuellement 
independants et uniformement distribues entre 0 et L 

On s'aper^oit qu'une simulation sequentielle est une operation S qui 
transforme un vecteur uniforme U= (Ui, U2, ....,Ursj ) en un vecteur structure 
Z = (Z], Z2,.... Zn) : 
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4 



Z = S(U) 



(1) 



Le probleme que represente la contrainte d'un vecteur Z a differents 
types de donnees peut etre resolu en contraignant les distributions 
conditionnelles Fc(Zi) (i = 1, N) et/ou le vecteur uniforme U= (Ui, U2, 



Des travaux recents sur Talgorithnie sequentiel ont ete concentres sur 
r amelioration de 1' estimation des distributions conditionnelles Fc(2j) par 
integration de donnees geologiques et de donnees sismiques. 

On peut citer par exemple le papier de Zhu, H. et al ; « Formatting and 
Integrating Soft Data : Stochastic Imaging via the Markov-Bayes Algorithm » 
in Soares, A., Ed. Geostatistics Troia 92, vol.1 : Kluwer Acad. Publ., 
Dordrecht, The Netheriands, p.1-12, 1993. 

Cependant, cette approche ne peut etre etendue a T integration de 
donnees non lineaires comme les pressions issues des tests de puits et les 
historiques de production, sauf a imposer une severe linearisation. En outre, 
comme toute combinaison de vecteurs uniformes U ne donne pas un vecteur 
uniforme, la methode de deformation graduelle d'un modele stochastique 
developpee dans le brevet deja cite, ne peut etre directement appliquee dans le 
cadre de la technique sequentielle rappelee ci-dessus. 

Aussi, la methode selon 1' invention permet de rendre compatible les 
deux approches i.e. d'etendre le formalisme developpe dans le brevet 
precedemment cite, a la deformation graduelle de realisations, generees par 
simulation sequentielle, d'un modele stochastique non necessairement 
gaussien. 
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Definition de la methode 

La methode permet de deformer graduellement une representation ou 
realisation, generee par simulation sequentielle, d'un modele stochastique non 
necessairement gaussien d'une grandeur physique z dans un milieu heterogene 
tel qu'une zone souterraine, afin de la contraindre a un ensemble de donnees 
collectees dans le milieu par des mesures et observations prealables, relatives a 
son etat ou sa structure. 

Elle est caracterisee en ce qu'elle comporte Tapplication d'un 
algorithme de deformation graduelle d'un modele stochastique a un vecteur 
gaussien (Y) ayant un nombre N de variables mutuellement independantes qui 
est rehe a un vecteur uniforme (U) a N variables uniformes mutuellement 
independantes par une fonction de repartition gaussienne (G), de fagon a 
definir une chaine de realisations u(t) du vecteur (U), et 1' utilisation de ces 
realisations u(t) pour generer des realisations z(t) de cette grandeur physique 
que Ton cale par rapport aux donnees (non lineaires). 

Suivant un premier mode de mise en ceuvre, on definit la chaine de 
realisations u(t) du vecteur uniforme (U) a partir d'une combinaison lineaire de 
realisations d'un vecteur gaussien (Y) dont les coefficients de combinaison 
sont tels que la sonmie de leurs carres est egale a 1' unite. 

Suivant un autre mode de mise en ceuvre, on realise une deformation 
graduelle d'un nombre n de parties du modele representatif du milieu 
heterogene tout en preservant la continuite entre ces n parties du modele, par 
une subdivision du vecteur uniforme (U) en un nombre n de sous-vecteurs 
mutuellement independants. 
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D'autres caracteristiques et avantages de la methode selon I'invention, 
apparaitront a la lecture de la description ci-apres d'un exemple non limitatif 
de realisation, en se referant aux dessins annexes ou : 

- la Fig.l montre la couche mediane d'une realisation d'un modele de facies 
servant de reference, generee par simulation sequentielle d'indicatrices ; 

- la Fig.2 montre la variation en fonction du temps de la pression obtenue 
dans un test de puits du modele de reference ; 

- les Fig. 3 A a 3E montrent respectivement cinq realisations initiales de la 
couche mediane d'une zone reservoir, contraintes seulement par les facies le 
long du puit ; 

- les Fig. 4A a 4E montrent respectivement pour ces cinq realisations, les 
courbes des pressions au puits dans le modele de reference comparees a 
celles correspondant aux modeles initiaux ; 

- les Fig.5A a 5E montrent respectivement cinq realisations de la couche 
mediane du modele de facies conditionnees aux facies le long du puits et 
calees par rapport a la courbe de pression obtenue par test de puits ; 

les Fig. 6A a 6E montrent respectivement pour les cinq reaUsations, les 
courbes des pressions au puits dans le modele de reference comparees a 
celles correspondant aux modeles cales ; 

- les Fig. 7 A a 7E montrent respectivement conmient les fonctions objectif 
correspondant respectivement a ces cinq exemples, varient avec le nombre 
d' iterations ; 
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- les Fig. 8 A a 8E montrent les deformations graduelles operees par un 
changement de coefficient d'anisotropie, sur un modele a trois facies genere 
par simulation sequentielle d'indicatrices; et 

- les Fig. 9A a 9E montrent les deformations graduelles locales d'un modele a 
5 trois facies, generees par simulation sequentielle d'indicatrices. 

Description detaillee de la methode 

On considere une zone d' etude que Ton subdivise par une grille a N 
mailles et Ton cherche a construire des realisations ou representations d'un 
modele stochastique d'une certaine grandeur physique z representant par 
10 exemple la permeabilite des formations dans la zone. Le modele recherche doit 
s'ajuster a des donnees obtenues par des mesures ou observations en un certain 
nombre de points, et notamment a des donnees non lineaires. 

Le calage d'un modele stochastique a des donnees non lineaires peut 
etre formulee comme un probleme d' optimisation. On designe par f*^^= (fi°^^ , 

15 f2°^^ , fs^^^ ... fp°^^), le vecteur des donnees non lineaires observees ou mesurees 

dans le milieu etudie (la zone reservoir), et par f = (fi , fa, , fp), le vecteur 

correspondant des reponses du modele stochastique de la permeabilite Z = (Zi, 
Z2,..., Zn). Le probleme de contraindre le modele stochastique Z par les 
observations consiste a generer une realisation z de Z qui reduit a une valeur 

20 assez faible une fonction objectif que Ton definit comme la somme des erreurs 
quadratiques ponderees des reponses du modele par rapport aux observations 
ou mesures dans la zone-reservoir, soit : 
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ou (Oi represente le poids attribue a la reponse fj. Les fonctiohs fj (i=l, 2, p) 
ainsi que la fonction objectif O sont des fonctions du vecteur Z. On se trouve 
devant un probleme d' optimisation de dimension N. 

Pour etendre le foraialisme developpe dans le brevet deja cite, a la 
5 deformation graduelle de realisations generees par simulation sequentielle non 
necessairement gaussienne, on va partir d'un vecteur gaussien a N variables Yi 
avec i = 1, 2, N, mutuellement independantes, de moyenne nulle et de 
variance egale a T unite, et on va definir N variables uniformes mutuellement 
independantes Ui, U2, U3,...Un par : 

10 Ui = G(Yi) V i = 1,2, ...,N 

oil G designe la fonction de repartition gaussienne standardisee. 

Ceci pose, on va appliquer Talgorithme de deformation graduelle 

developpe dans un cadre gaussien, au vecteur gaussien Y = (Yi, Y2, , Yn), 

pour constmire une chaine continue de realisations du vecteur uniforme U = 

15 (Ui, U2, Un). Etant donne deux realisations independantes yaCt yb de Y, on 

definit la chaine de realisations u(t) du vecteur U que Ton obtient selon la 
relation : 

u(t) = Giy^ cost-^y, sin/) (2) 

Pour chaque t, u(t) est une realisation du vecteur U. Ensuite, par 

20 echantillonnage de la distribution conditioimelle F^(z.) (i=l, 2, N) en 

utilisant les elements du vecteur u(t), on obtient un vecteur z(t) qui est, pour 
chaque t, une realisation du vecteur aleatoire Z. En consequence, on pent 
ajuster le parametre t comme dans le cas gaussien de fa9on a caler z(t) a des 
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donnees non lineaires. Cette procedure est iteree jusqu'a atteindre un calage 
satisfaisant. 

Calage d'un modele de fades a des donnees de pression obtenues 
par des tests de puits 

5 Pour illustrer T application de la methode d' optimisation stochastique ci- 

dessus definie, on va chercher a caler un modele stochastique de reservoir a 
des donnees de pressions obtenues par des tests de puits. La construction du 
modele de reservoir est derivee d'une formation petroliere reelle qui comporte 
trois types de facies : deux facies reservoir de bonne qualite (facies 1 et 2) et 

10 un type de facies reservoir de tres mauvaise qualite (facies 3). La table 1 
definit les proprietes petrophysiques des trois facies : 







K,{md) 


K.{md) 




c,{\0-^bar-') 


facies 1 


10. 


10. 


10. 


17 


2.1857 


facies 2 


1. 


1. 


1. 


14 


2.0003 


facies 3 


0.1 


0.1 


0.001 


9 


1.8148 



Pour representer la distribution specifique de facies de la formation 
petrolifere, on genere d'abord une realisation binaire pour representer le facies 
15 3 et son complement. Ensuite, dans la partie complementaire du facies 3, on 
genere une autre realisation binaire independante de la premiere pour 
representer les facies 1 et 2. La formation est discretisee par un maillage 
regulier de 60x59x15 blocs de taille 15mxl5mxL5m. Un modele de 
variogramme exponentiel est utilise pour estimer les distributions 
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conditionnelles. La direction d' anisotropic principale est diagonale par rapport 
au maillage. Les portees du variogramme du facies 3 selon les trois directions 
d'anisotropie sont respectivement de 300m, 80m et 3m. Les portees du 
variogramme des facies 1 et 2 selon les trois directions d'anisotropie sont 
5 respectivement de 150m, 40m et L5m. Les proportions des facies 1, 2, 3 sont 
respectivement de 6%, 16% et 78%. 

On a realise un test de puits au moyen d'un simulateur de test de puits a 
differences finies, de la maniere decrite par : 

Blanc, G. et al : « Building Geostatistical Models Constrained by Dynamic 
10 Data - A Posteriori Constraints » in SPE 35478, Proc. NPF/SPE European 
3D Reservoir Modelling Conference, Stavenger, Norway, 1996. 

On voit sur la Fig.l la couche mediane d'une realisation qui sert de 
modele de reference pour cette validation. La section du puits qui a ete 
perforce traverse horizontalement la couche mediane du modele de reservoir le 
15 long de Taxe x. Le diametre du puits est de 7.85cm, la capacite du puits est 
nulle et les facteurs parietaux (skin) des facies 1, 2 et 3 sont respectivement 0, 
3 et 50. Le test de puits synthetique dure 240 jours avec un debit constant de 5 
m3/jour de fagon a investiguer presque la totalite du champ petrolifere. La 
Fig.2 montre la variation de pression en fonction du temps. 

20 Le but etait de construire des realisations du modele de facies contraint 

par les facies rencontre le long du puits et par la courbe de pression obtenue 
lors du test de puits. La fonction objective est definie conune la somme des 
differences quadratiques entre les reponses en pression du modele de reference 
et les reponses en pression de la realisation. Du fait que le comportement 

25 dynamique du modele de reservoir est principalement controle par le contraste 
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entre les facies reservoir de bonne et de mauvaise qualites, la realisation 
binaire utilisee pour generer les facies 1 et 2 a ete fixee d'abord et seule la 
realisation binaire utilisee pour generer le facies 3 a ete deformee pour le 
calage des donnees de pression. 

5 Les reponses en pression resultant des tests de puits pour les cinq 

realisations des Fig.3A a 3E, sont differentes de celle du modele de reference, 
comme le montrent les Fig. 4A a 4E. Partant respectivement de ces 5 
realisations independantes, en utilisant la methode de calage iteratif ci-dessus, 
on obtient apres quelques iterations, cinq realisations calees (Fig. 5 A a 5E) 
10 pour lesquelles les courbes de pression correspondantes sont tout a fait en 
accord avec celles du modele de reference, conune le montrent les Fig.6A a 
6E. 

Deformation graduelle par rapport aux parametres structuraux 

Dans de nombreux cas, on ne dispose pas de doimees suffisantes pour 
15 en deduire les parametres structuraux du modele stochastique : moyenne, 
variance, fonction de covariance, etc. Ces parametres structuraux sont sou vent 
donnes en termes d'intervalles ou de distributions a priori. Si leurs valeurs sont 
fausses, il est vain de rechercher une realisation calee sur des donnees non 
lineaires. Aussi, il est essentiel pour les applications de pouvoir proceder a une 
20 deformation graduelle d'une realisation avec modification simultanee de 
nombres aleatoires et des parametres structuraux. L'algorithme de simulation 
sequentielle definie par T equation (1) foumit cette possibilite de changer 
simultanement ou separement Toperateur structurel S et le vecteur uni forme U. 
Les Fig.8A-8E montrent les deformations graduelles que Ton obtient pour une 
25 realisation fixe du vecteur uniforme U quand on change le coefficient 
d' anisotropic. 



wo 01/02876 




PCT/FROO/01853 



Deformation graduelle locale ou regionalisee 

Quand les observations sont reparties sur des zones differentes d'une 
formation etudiee, un calage qui utiliserait une deformation globale risquerait 
d'etre inefficace car Tamelioration de I'accord que Ton obtiendrait sur une 
zone pourrait le deteriorer sur une autre zone. II est done bien preferable 
d'appliquer des deformations graduelles zone par zone. On considere une 

1 2 

subdivision du vecteur U en un certain nombre n de sous-vecteurs U , U 
U" mutuellement independants ce qui rend possible de realiser 
individuellement leur deformation graduelle. En appliquant separement a 
chacun des sous-vecteurs U\ U^,...., U", I'algorithme de deformation 
graduelle, on obtient une fonction de dimension n, du vecteur uniforme U : 

G(rJ cosf, -f sin^t) 
G(y/ cos r 2 + y/sin t^) 

G(y; cos + sin fj_ 

ou yj et y/ pour tout i =1, 2, n, sont des sous-vecteurs gaussiens 
independants. Pour un ensemble de realisations de yj et , on resout un 
probleme d' optimisation de n parametres tj, t2, tn pour obtenir une 
realisation qui ameliore, sinon maintient, le calage des donnees. Cette 
procedure pent etre iteree jusqu'a atteindre un calage satisf aisant. 

En realisant des deformations locales graduelles, on ameliore de fa9on 
significative la vitesse de calage dans tous les cas oil les mesures ou 
observations sont reparties dans differentes zones du milieu. 
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On voit bien I'effet de cette deformation locale graduelle sur le modele 
a trois facies, des Fig. 9 A a 9E ou seule la partie inferieure gauche delimitee 
est af fee tee. 

La methode selon T invention pent etre facilement generalisee a la 
5 deformation graduelle d'une representation ou realisation d'un modele 
stochastique quelconque, car generer une realisation d'un tel modele 
stochastique se ramene toujours a la generation des nombres uniformes. 
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REVENDICATIONS 

1) Methode pour deformer graduellement une representation ou 
realisation, generee par simulation sequentielle, d'un modele stochastique non 

5 necessairement gaussien d'une grandeur physique z dans un milieu heterogene 
tel qu'une zone souterraine, afin de la contraindre a un ensemble de donnees 
coUectees dans le milieu par des mesures et observations prealables, relatives a 
son etat ou sa structure, caracterisee en ce qu'elle comporte T application d'un 
algorithme de deformation graduelle d'un modele stochastique a un vecteur 

10 gaussien (Y) a N variables mutuellement independantes qui est relie a un 
vecteur uniforme U a N variables uniformes mutuellement independantes par 
une fonction de repartition gaussienne (G), de fagon a construire une chaine de 
realisations u(t) du vecteur U, et Tutilisation de ces realisations u(t) pour 
generer des realisations z(t) de cette grandeur physique, que Ton cale par 

15 rapport aux donnees. 

2) Methode selon la revendication 1, caracterisee en ce que Ton definit 
une chaine de realisations u(t) du vecteur (U) a partir d'une combinaison 
lineaire de realisations du vecteur gaussien (Y) dont les coefficients de 
combinaison sont tels que la somme de leurs carres est egale a T unite. 

20 3) Methode selon la revendication 1 ou 2, caracterisee en ce qu'elle 

comporte une deformation graduelle du modele representatif du milieu 
heterogene, simultanement par rapport aux parametres structuraux et aux 
nombres aleatoires. 

4) Methode selon la revendication 1 ou 2, caracterisee en ce qu'elle 
25 comporte une deformation graduelle separee d'un nombre n de parties du 
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modele representatif du milieu heterogene tout en preservant la continuite 
entre ces n parties du modele, par une subdivision du vecteur unifomie en n 
sous-vecteurs mutuellement independants. 
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